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Abstract. We develop Hamilton-Jacobi theory for Chaplygin systems, a certain class of non- 
holonomic mechanical systems with symmetries, using a technique called Hamiltonization, which 
^ transforms nonholonomic systems into Hamiltonian systems. We give a geometric account of the 

^ Hamiltonization, identify necessary and sufficient conditions for Hamiltonization, and apply the 

^— ■») conventional Hamilton-Jacobi theory to the Hamiltonized systems. We show, under a certain suffi- 

cient condition for Hamiltonization, that the solutions to the Hamilton-Jacobi equation associated 
, ^ with the Hamiltonized system also solve the nonholonomic Hamilton-Jacobi equation associated 

with the original Chaplygin system. The results are illustrated through several examples. 

o 

^ 1. Introduction 

' I ' 1.1. Background and Motivation. In 1911 S.A. Chaplygin published a paper (re-published in 

pi I English in [8j) introducing his theory of the "reducing multiplier" into the study of nonholonomi- 

cally constrained mechanical systems. In his paper, Chaplygin showed that a two degree of freedom 
nonholonomic system possessing an invariant measure became Hamiltonian after a suitable repa- 
rameterization of time, a process we would like to refer to as Chaplygin Hamiltonization. Since 
I ^ I then, Chaplygin's result has generated considerable interest and been extended [iTl [HI [151 EHl EH] 

to more general settings. 

However, a second contribution contained in Chaplygin's paper has been left undeveloped. In 
Section 5 of his paper, Chaplygin integrates the nonholonomic system now known as the Chaplygin 
\^ Sleigh [3j by using the Hamilton-Jacobi equation for the Hamiltonized system. The aim of this 

paper is to develop this idea further to establish a link with the nonholonomic Hamilton-Jacobi 
equation in Iglesias-Ponte et al. [16] and Ohsawa and Bloch [25j. 

Specifically, we first employ the technique called Chaplygin Hamiltonization to transform Chap- 
I lygin systems into Hamiltonian systems, and then apply the conventional Hamilton-Jacobi theory 

T-H to the resulting Hamiltonian systems to obtain what we would like to call the Chaplygin Hamilton- 

Jacobi equation. This is an indirect approach towards Hamilton-Jacobi theory for nonholonomic 
• '-j systems, compared to the direct approach of extending Hamilton-Jacobi theory to nonholonomic 

systems, as in Iglesias-Ponte et al. |16j . de Leon et al. [H], Ohsawa and Bloch [23], and Cariiiena 
c5 et al. |7]. 

1.2. Direct vs. Indirect Approaches. The indirect approach to nonholonomic Hamilton-Jacobi 
theory via Chaplygin Hamiltonization has both advantages and disadvantages. The main advantage 
is that we have a conventional Hamilton-Jacobi equation and thus the separation of variables 
argument applies in a rather straightforward manner compared to the direct approach in Ohsawa 
and Bloch [25]. A disadvantage is that the Chaplygin Hamiltonization works only for certain 
nonholonomic systems; even if it does, the relationship between the Hamilton-Jacobi equation 
and the original nonholonomic system is not transparent, since one has to inverse-transform the 
information in the Hamiltonized systems. Nevertheless, Hamiltonization is known to be a powerful 
technique for integration of nonholonomic systems [H [SJ [TT| [T2[ [H] , and hence it is interesting to 
establish a connection with the direct approach. 
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Let us briefly summarize the differences between two approaches. Recall from Ohsawa and Bloch 
[25j that the nonholonomic Hamilton- J acohi equation is an equation for a one-form 7 on the original 
configuration manifold Q: 

Hoj = E, (1.1) 

along with the condition that 7, seen as a map from Q to T*Q, takes values in the constrained 
momentum space 7W c T*Q (see Eq. (2.4) below), i.e., 7 : Q — -M, and also that 

dllvxv = 0, i.e., d'ylv^w) = for any v,w gV, (1-2) 

where T> C TQ is the distribution defined by nonholonomic constraints, and H : T*Q ^ M the 
Hamiltonian. 

On the other hand, the Chaplygin Hamiltonization first reduces the system by identifying it 
as a so-called Chaplygin system with a symmetry group G, and then Hamiltonizes the system on 
the cotangent bundle T*{Q/G) of the reduced configuration space Q/G. The resulting system is 
a (strictly) Hamiltonian system on T*(Q/G) with another Hamiltonian Hq : T*(Q/G) — )• M; so 
we may apply the conventional Hamilton- Jacobi theory to the Hamiltonized system to obtain the 
Chaplygin Hamilton- J acohi equation 

HcodW = E, 

which is a partial differential equation for a function W : Q/G — )• M. Therefore, the difference 
lies not only in the forms of the equations (the former involves the one-form 7, which is not even 
closed, whereas the latter invokes the exact one- form dW), but also in the spaces on which the 
equations are defined. Furthermore, the Chaplygin Hamilton-Jacobi equation corresponds to the 
Hamiltonized dynamics and is related to the original nonholonomic one in a rather indirect way. 
Therefore, on the surface, there does not seem to be an apparent relationship between the two 
approaches. 

1.3. Main Results. The main goal of this paper is to establish a link between the two distinct 
approaches towards Hamilton-Jacobi theory for nonholonomic systems. To that end, we first 
formulate the Chaplygin Hamiltonization in an intrinsic manner to elucidate the geometry involved 
in the Hamiltonization. This gives a slight generalization of the Chaplygin Hamiltonization by 
Fedorov and Jovanovic fTT] and also an intrinsic account of the necessary and sufficient condition 
for Hamiltonizing a Chaplygin system presented in [IHj. These results are also related to the 
existence of an invariant measure in nonholonomic systems (see, e.g., Kozlov [22], Zenkov and 
Bloch f30j, and Fedorov and Jovanovic [llj). 

We also identify a sufficient condition for the Chaplygin Hamiltonization, which turns out to be 
identical to one of those for another kind of Hamiltonization (which renders the systems "conformal 
symplectic" [l5]) obtained by Stanchenko [28] and Cantrijn et al. [6]. We then give an explicit 
formula that transforms the solutions of the Chaplygin Hamilton-Jacobi equation into those of the 
nonholonomic Hamilton-Jacobi equation (see Fig. [T]). Interestingly, it turns out that the sufficient 



Chaplygin System 



Hamiltonization 



Nonholonomic 
H- J Theory [l6l [25] 



H = E, c?7|x)xD = 



Hamiltonized Chaplygin System 



H-J Theory 



HcodW = E 



Figure 1. Relationship between the nonholonomic H-J equation applied to a Chap- 
lygin system and the H-J equation applied to the Hamiltonized Chaplygin system 
Explicit formulas for the correspondence W i— )• 7 are given in Theorems 4.1 and 7.1 
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condition plays an important role here as well. We also present an extension of these results to 
a class of systems that are Hamiltonizable after reduction by two stages, following the idea of 
Hochgerner and Garcfa-Naranjo [15j. We illustrate, through several examples, that the Chaplygin 
Hamilton- Jacobi equation may be solved by separation of variables, and that the solutions are 
identical to those obtained by Ohsawa and Bloch [25] after the transformation mentioned above. 



1.4. Outline. We begin with an overview of nonholonomic mechanical systems in Section 2.1 



specializing to Chaplygin systems in Section 2.2 After discussing the relationship between the 



Hamiltonizability of a nonholonomic system and the existence of an invariant measure for it in Sec- 



tion 3.1 , we derive necessary and sufficient conditions for a Chaplygin system to be Hamiltonizable 



in an intrinsic manner in Section 3.2 This result then leads to the development of Hamilton- Jacobi 



theory for Hamiltonizable Chaplygin systems in Section 4.2, Specifically, we relate the Chaplygin 
Hamilton- Jacobi equation for the Hamiltonized system with the nonholonomic Hamilton-Jacobi 
equation for the original system. A couple of examples are presented in Section [5] to illustrate 
the theoretical results. In Section [6j we introduce a further reduction of the reduced Chaplygin 
systems under certain conditions; the second reduction is employed to Hamiltonize those systems 
that are not Hamiltonizable after the first reduction. Then, in Section [7| we relate the Chaplygin 
Hamilton-Jacobi equation for such systems with the nonholonomic Hamilton-Jacobi equation. We 
then illustrate the theory in the Snakeboard example. 

2. Chaplygin Systems 

2.1. Hamiltonian Formulation of Nonholonomic Mechanics. Consider a nonholonomic sys- 
tem on an n-dimensional configuration manifold Q with a constraint distribution T) C TQ defined 
by the constraint one-forms {w^}^^ as 

V = {v e TQ \ u}'{v) = Q, s = I, . . . ,m} 

and also with the Lagrangian L : TQ — t- M of the form 

Livg) = ^gg{v„vg)-V{q) (2.1) 



with the kinetic energy metric g defined on Q. Define the Legendre transform ¥L : TQ T*Q hy 



{¥L{Vq),Wq) = ^L{Vq + eWq) 



e=0 

where the last equality defines g^ : TQ T*Q; hence we have ¥L = g^ . Also define the Hamiltonian 
H : T*Q ^ M by 

H{Pq) ■■= {Pq,Vq) - L{Vq), 

where Vq = (FL)^^(pg) on the right-hand side. Then, Hamilton's equations for nonholonomic 
systems are written as follows: 

ix^ = dH - Xs-k*qUj' , (2.2) 

along with 

Tttq{X)£V or u'iTiTQiX)) = for s = l,...,m, (2.3) 

where ttq : T*Q ^ Q is the cotangent bundle projection. Introducing the constrained momentum 
space 

M ■.= ¥L{V) CT*Q, (2.4) 
the above constraints may be replaced hy p £ Ai. 
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2.2. Chaplygin Systems. 

Definition 2.1 (Chaplygin Systems). A nonholonomic system with Hamiltonian H and distribu- 
tion D is called a Chaplygin system if there exists a Lie group G and a free and proper group action 
of it on Q, i.e., $ : G x Q — > Q or : Q — )• Q for any h G G, such that 

(i) the Hamiltonian H and the distribution D are invariant under the G-action; 

(ii) for each q G Q, the tangent space TqQ is the direct sum of the constraint distribution and the 
tangent space to the orbit of the group action, i.e., 

TqQ =Vg® TqOg, 

where Oq is the orbit through q of the G-action on Q, i.e., 

Oq := {Mq) GQ\hGG}. 
This setup gives rise to the principal bundle 

TT-.Q^ Q/G =: Q 

and the connection 

A-.TQ^g, 

with being the Lie algebra of G such that ker A = V. So the above decomposition may be written 

as 

TqQ = kei Aq e ker Tgvr. 

Furthermore, for any q G Q and q := Tr{q) S Q, the map TIj7r|x), : Vq — )• TqQ is a linear isomorphism, 
and hence we have the horizontal lift 

hl^ : TqQ ^ Vq- Vq ^ (T,^| 2, J " ^ . 

We will occasionally use the following shorthand notation for horizontal lifts: 

:= hl^iVq). 



Therefore, any vector Wq € TqQ can be decomposed into the horizontal and vertical parts as follows: 

Wq = hor(Wg) + ver(Wg), (2.5a) 

with 

hor(M^,) = hl^{Wq), VeiiWq) = {Aq{Wq))Q{q), (2.5b) 

where Wq := TqTT{Wq) and G ^{Q) is the infinitesimal generator of G g. 
We may then define the reduced Lagrangian 

L:=Lohf, (2.6a) 

or more explicitly, 

L:TQ^R; Vq ^ ^gq{vq, Vq) - V{q), (2.6b) 
where g is the metric on the reduced space Q induced by g as follows: 

9q{Vq,Wg) := gg{h\g {Vq),hlg (Wg)} = gg{Vq,Wq), (2.7) 

and the reduced potential y : Q — ?• M is defined such that V = V o tt. 

This geometric structure is carried over to the Hamiltonian side (see Ehlers et al. [IQJ). Specifi- 
cally, we define the horizontal lift hl^ : T^Q ^ Mg by 

hlf := ¥Lg o h^ o(FL)ri = g^ o hl^ o{-g%\ (2.8) 
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or the diagram below commutes. 



I <3 



T-nQ 



We will use the shorthand notation 



hlf (a,) 
for any aq G T|Q. 

We also define the reduced Hamiltonian H : T*Q 



by 



(2.9) 

It is easy to check that this definition coincides with the following one by using the reduced La- 
grangian L: 

Hijpq) ■■= {Vq,Vq) - L{vq), 

with Vq = {¥L)q^{pq). 

Performing the nonholonomic reduction of Koiller [18] (see also Bates and Sniatycki p] , Ehlers 
et al. [10], and Hochgerner and Garci'a-Naranjo |15j). we obtain the reduced Hamilton's equations 
for Chaplygin systems defined by 







nh 



dH 



(2.10) 

with the almost symplectic form 

Q^'^:=^-'~, (2.11) 

where X is a vector field on T*Q and Cl is the standard symplectic form on T*Q; the two- form H 
on T*Q is defined as follows: For any aq G T^Q and yan-,Zan G Ta,T*Q, let Yq := ^71(5(3^0 ) and 



Zq := TiTQ^Zaq) where ttq : T*Q ^ Q is the cotangent bundle projection, and then set 

H„,(3^„„2.,) := (Johlf (a5),^,(hl^(y^),hl^(Z,-))) 

= {j{a^^),BqiY^,Z^)), (2.12) 

where J : r*Q — )• g* is the momentum map corresponding to the G-action, and B is the curvature 
two-form of the connection A. This is well-defined, since the Ad*-equivariance of the momentum 
map J and the Ad-equivariance of the curvature B cancel each other [19j: Writing hq := <I>/j(g), we 



have, using Lemma A.l and the G-equivariance of the momentum map J and the curvature B (see. 



e.g., Marsden et al. [211 Corollary 2.1.11] for the latter), 

;j«),fi,,(y4,z,\)) = (j(r*$,-i(a^)) ,c^;;fi,(y^^z, 

JAd^,J(a^),Ad,i3,(y,^Z^> 
[j{a^^),Bq{Y^\Z^ 



3. Chaplygin Hamiltonization of Nonholonomic Systems 
This section discusses the so-called Chaplygin Hamiltonization of the reduced dynamics defined 



by Eq. (2.10). The results here are mostly a summary of some of the key results of Stanchenko 

However, our 



[28], Cantrijn et al. |6], Fedorov and Jovanovic [llj, and Fernandez et al. [E 
exposition is slightly different from them, and also touches on those aspects that are not found 
in the above papers. Furthermore, our intrinsic account of the Hamiltonization provides us with 
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a better understanding of the geometry involved in it, and then leads us to our main results on 
nonholonomic Hamilton-Jacobi theory in Sections |4] and [7j 



3.1. Hamiltonization and Existence of Invariant Measure. We first discuss the relationship 
between Hamiltonization and existence of an invariant measure for nonholonomic systems. The 



next subsection will show how to Hamiltonize the reduced system, Eq. (2.10), explicitly. 

Let / : T*(5 — 5- M be a smooth nowhere- vanishing function that is constant on each fiber, i.e., 
f{otq) = fifiq) for any Og, f3q G T^Q- Therefore, we can write, with a slight abuse of notation, 
f{oiq) = f{q); so / may be seen as a function on Q as well. 

Remark 3.1. The above definition of the function / is essentially the same as that of Chaplygin 
[8], where / is defined as a function on Q. However, in the present work, it is more convenient to 
formally define / function on T*Q. 

Remark 3.2. In the discussion to follow, we derive certain conditions on the function / in order 



to Hamiltonize the system given by Eq. (2.10). It sometimes turns out that such / is nowhere- 



vanishing only on an open subset U in Q. In such cases, we redefine Q := U. 
Now, consider the vector field 

1 



X/f = jXeX{T*Q), 



X / f 

and let <I>j : T*Q — )• T*Q be the flow defined by the corresponding vector field, i.e., for any 



iX/f){ag 



t=o 



Xiag). 



Furthermore, define a map : T*Q ^ T*Q hy 

^'Z : a I— )• fa, 

which is clearly a diffeomorphism with the inverse "^^^ = ^'i/j : T*Q T*Q; a ^ a/f, and define 
$f c . T*Q _^ T*Q by 



or the diagram below commutes. 



x/f 



T*Q 
*7'='fi// 

T*Q 



T*Q 



T*Q 



a// I -$f/^(a//) 



a \- 



Y 



a) 



(3.1) 
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Then, we have the vector field Xq € X{T*Q) corresponding to the flow which is the pull-back 
of X/f by = ^yf. For any a,- G T*Q, 

d 



XcK-) := ^'^f^K; 



t=o 



t=Q 



= T^f{X/my\a,)) 
= i^f'nX/f)ia,) 

= ^*/^(X//)K). (3.2) 

In particular, the third line in the above equation shows that X/f and Xq are ^'/-related: 

r% o (X/f) = Xc o (3.3) 

Now, we relate relate the (possible) symplecticity of the vector field Xq with the existence of an 
invariant measure for the reduced system, Eq. (2.10): 



Theorem 3.3. If Xq G X{T*Q) is symplectic, i.e., £xq^ = 0? then the reduced system, Eq. (2.10), 
has the invariant measure /""^A, where n := dimQ and A is the Liouville volume form 

A := ^ '— p A • - A 9. = • • • A dg" A dpi A • • • A dpn- 

n 

In other words, we have 

i:^(r-iA) = o. 

This theorem is a slight generalization of the following: 
Corollary 3.4 (Fedorov and Jovanovic [H]). If Xq G X{T*Q) is Hamiltonian, i.e., 

ix^^ = dHc 

for some Hq : T*Q — )• R, then the reduced nonholonomic dynamics, Eq. (2.10), has the invariant 
measure f^~^A. 

Proof. Follows easily from Cartan's formula: 

£^^n = diix^d) + ix^dn = ddHc = 0. □ 

We state a couple of lemmas before proving Theorem |3.3[ 

Lemma 3.5. Let f : T*Q — t- M 6e a smooth function that is constant on each fiber, i.e., f{ceq) = 
f{/3g) for any aq,(3g £TiQ. Then, 

{^}n) A • • • A = f"nA-- - An. 

' ^ ' ^ n ' 

n 

Proof. Let O be the symplectic one-form on T*Q, i.e., 0. = —dQ. Let us first calculate ^JQ: We 
have, for any a & T*Q and v G TaT*Q, 

(*;e),(t;) = e^^(„)(T%(i;)) 

= {^f{a),TTTQoT^f{v)) 
= {fa,T{7TQO^^>f){v)) 
= f {a,T7rQ{v)) 

= feaiv), 
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where we used the fact that ^/ is fiber-preserving, i.e., ttq o "^j = ttq. Hence we have ^'^0 = /9, 
and thus 

= -dim 

= -df A e - fde 

= fn-dfAe. (3.4) 

Therefore, using the fact that a A P = (3 A a for any two- forms a and /3, we have 

(^p) A • • • A = rCiA-'-An 

+ Cf] (_i)fc/"-fc pA- - - aO a (d/ a g) a • • • a (d/ a e) . 

Let us show that the second term vanishes. Since / is constant on fibers, we have 



Therefore, 



df AQ=pb^dq'' Adq' 



and thus df A Q does not contain any term with dpaS. On the other hand, fl A ■ — A i\ contains 

n—k 

only n — A; of dpaS. Therefore, the 2n-form 

A • • • A A (4f A G) A • • • A (d/ A 9) 
' ^ ^ ' 

n—k 1^ 

contains only n — k of dp^s, and thus n+k of dg'^'s, which implies that this 2n-form must vanish. □ 

Definition 3.6. Let M be an n-dimensional orientable manifold, and /u be a volume form, i.e., a 
nowhere-vanishing n-form. Then, the divergence div^(X) of a vector field X on M relative to fi is 
defined by 

fxfJ- = div^{X) fi. (3.5) 
Therefore, the flow of X is volume- preserving if and only if div^{X) = 0. 

Lemma 3.7. Let M be an orientable differentiable manifold with a volume form fi, X a vector 
field on M , and f a nowhere-vanishing smooth function on M . Then, the following identity holds: 

div^(/X) = /divj^(X). (3.6) 

Proof. We have the identities [see, e.g.,[l, Proposition 2.5.23 on p. 130] 

divy^(X) = div^(X) + jX[fl div^ifX) = fdiv^{X) + X[f]. 

Multiplying the first by / and taking the difference of both sides, we have the desired identity. □ 



Proof of Theorem 3.3. As shown in Eq. (3.3), the vector fields X / f and Xq are ^'/-related. There- 
fore, 
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since Xq is assumed to be symplectic; thus 



However, by Lemma 3.5 we have 



and hence £jf/j(/"A) = 0; this imphes divjfi^(X//) = 0. Then, the above lemma gives 

div^.-i^W = /div/.AW/) = 0, 

which imphes £x{P~^^ = 0. □ 

3.2. The Chaplygin Hamiltonization. Here we discuss the so-cahed Chaplygin Hamiltonization 
of the reduced system, Eq. (2.10). Let us first find the equation satisfied by the vector field Xq 
defined in Eq. (3.2). 

Lemma 3.8. The vector field Xq ^ X(T*Q) satisfies the following equation: 



where Hc:T*Q 



is defined by 



1 



n + j{dfAe-fE) 



He ■= H o ^i/f. 



dHc, 



(3.7) 



(3.8) 



Proof. As shown in Eq. (3.3), the vector fields X/ f and Xc are ^'j-related. Therefore, ^'rij^^a 



have 



'^a for any differential form a [see, e.g.,|Tl Proposition 2.4.14]; in particular, for a = fi, we 



Using Eqs. (3.4) and (2.10) on the right-hand side, we have 

ix/f^}n = i^/f{fn-dfAQ) 

= ix^-ix/f{df AO) 

= dH + ixE-ix/f{df A&) 

= dH-ix/f{dfAe-fE). 

Therefore, 

'^}ixc^ + ix/f{dfAe-fE)=dH, 
and then applying to both sides gives 



ix^n + ^l^fix/f{df AO - fE) = dHc. 



Since the vector fields Xq and X/f are ^'i/j-related, we have '^*^^j:ix/fO: = ix^'^l/jCt for any 
differential form a; hence 



1// ''X/f 



A e - / H) = ixan/f{df A e - / H) 



'Xc 



di^l/ff)Ai^*yf@)-^*yfifE 



1// 



^Xc[#A(G//)-S] 



'Xc 



f 



{dfAQ-fE) 
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0/ / as in the proof of Lemma 



3.5 



where ^i/ff = f since / is constant on each fiber; ^*/j0 
^yj'E = H// follows from the following calculation: From the definition of S in Eq. (2.12), we have 

( J o hlf (ag/f) , B, {hlf (Yg) , hl^ (Z^)) > 
^ (johl^(a5),^,(hl^(y,),hl^(Z,))> 



fiQ) 
1 

m 



tqiyaql ^Qq) 



where, in the second line, we defined Yq, Zq G TqQ as 

Yq := Tttq o T^i/f{ya,) = T(7rQ o ^i/f){ya,) = TTrgiy^^), 
and Zq in the same way, which coincide the ones introduced earlier when defining H; the third line 
follows from the linearity of hl"^ and also of J in the fiber variables. □ 

Proposition 3.9 (Necessary and Sufficient Condition for Hamiltonization) . The vector field Xq G 
X{T*Q) satisfies Hamilton's equations 

ix^^ = dHc (3.9) 
if and only if the one-form i^^ (df f\Q — f vanishes. 

Proof. Follows immediately from Lemma |3.8[ □ 

Remark 3.10. Locally, the above necessary and sufficient condition is precisely Eq. (2.17) in Fer- 
nandez et al. |l4j . 

Definition 3.11. The process of finding an / satisfying the above condition is called Chaplygin 



Hamiltonization, or just Hamiltonization for short; the resulting Hamiltonian system, Eq. (3.9), is 
called the Hamiltonized system; we would like to call Hq a Chaplygin Hamiltonian. 



Now, combining Proposition 3.9 with Theorem 3.3 or Corollary 3.4, we have 



Corollary 3.12. Suppose there exists a nowhere-vanishing fiber-wise constant function f : T*Q — t- 
M such that ix^{df A — / H) vanishes. Then, the 2n-form f^~^A is an invariant measure of the 
reduced system, Eq. (12.10). 



We now state the main result of this section. The following theorem will be used in the next 
section in relation to the nonholonomic Hamilton-Jacobi theory: 

Theorem 3.13 (A Sufficient Condition for Hamiltonization). Suppose there exists a nowhere- 
vanishing fiber- wise constant function f : T*Q ^M. that satisfies the equation 

dfA@ = fE. (3.10) 

Then, the vector field Xq € X(T*Q) satisfies the following Hamilton's equations: 



ix^i} = dHc, 



(3.11) 



and, as a result, the reduced nonholonomic dynamics Eq. (2.10) has the invariant measure /" ^A. 
Proof. Straightforward from Lemma 3.8 and Corollary 3.4, □ 



Remark 3.14. Locally, the sufficient condition (3.10) becomes condition (2.22) in Fernandez et al. 
T4i. 



Remark 3.15. As shown by Stanchenko [2^ (see also Cantrijn et al. |6J), Eq. (3.10) is also a sufficient 
condition for the two-form Clf := /(O — H) to be closed, so that Eq. (2.10) becomes 

ix/fClf = dS, 

and so the dynamics of X// is Hamiltonian with the non-standard symplectic form j . 
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4. NONHOLONOMIC HAMILTON-JACOBI THEORY VIA ChAPLYGIN HAMILTONIZATION 



4.1. The Chaplygin Hamilton— Jacobi Equation. Since the Hamiltonized system, Eq. (3.11) 



is a canonical Hamiltonian system on T*Q, we may apply the conventional Hamilton- Jacobi theory 
(see, e.g., Abraham and Marsden [1, Chapter 5]) to the system and obtain the (time-independent) 
Hamilton- Jacobi equation: 

HcodW = E, (4.1) 
with an unknown function W : Q ^ M. and a constant E (the total energy). We would like to call 



Eq. (4.1) the Chaplygin Hamilton-Jacobi equation. 



Now that we have two Hamilton-Jacobi equations for Chaplygin systems, i.e., the nonholonomic 



Hamilton- Jacobi equation ( 1.1 ) and the Chaplygin Hamilton-Jacobi equation (4.1 ), a natural ques- 



tion to ask is: What is the relationship between the two? 
4.2. Relationship between the Chaplygin H J and Nonholonomic H J Equations. In 



relating the Chaplygin Hamilton-Jacobi equation (4.1) to the nonholonomic Hamilton-Jacobi equa- 



tion (1.1), a natural starting point is to look into the relationship between the Chaplygin Hamil- 



tonian He and the original Hamiltonian H (recall from Eqs. (2.9) and (3.8) that they are related 



through the Hamiltonian H); the upper half of the following commutative diagram shows their 
relationship. 




T*Q 



(4.2) 



dW 



Now, suppose that a function 1^ : Q — )• M satisfies the Chaplygin Hamilton-Jacobi equation (4.1). 
This means that the one-form dW, seen as a map from Q to T*Q, satisfies Hq ° dW{q) = E for 
any q G Q with some constant E; equivalently, Hq ° dW o '7T{q) = E for any q G Q. The lower half 
of the above diagram (4.2) incorporates this view, and also leads us to the following: 



Theorem 4.1. Suppose that there exists a nowhere-vanishing fiber-wise constant function f 
T*Q ^ M that satisfies Eq. (3.10), and hence by Theorem 
( 3.11[ ) for the vector field Xc- Let W : Q - 
equation (4.1), and define 7 : Q — )■ 7W by 



3.13, we have Hamilton's equations 



be a solution of the Chaplygin Hamilton-Jacobi 



-f{q) := hl^ 0^1/^ odWo Tr{q) 



hi 



M 



1 



m 



dW{q) 



(4.3) 



where q := 7r{q). Then 7 satisfies the nonholonomic Hamilton-Jacobi equation (1.1) as well as the 
condition Eq. (|1.2|). 



Remark 4.2. Notice that Theorem 4.1 relates a solution of the Chaplygin Hamilton-Jacobi equation. 



which is for the reduced dynamics defined by Eq. (3.11 ), with that of the nonholonomic Hamilton 



Jacobi equation for the full dynamics defined by Eq. (2.2). Therefore, the theorem provides a 



method to integrate the full dynamics by solving a Hamilton-Jacobi equation for the reduced 
dynamics. 



Proof. That the one-form 7 defined by Eq. (4.3) satisfies the nonholonomic Hamilton-Jacobi equa- 



tion (1.1) follows from the diagram (4.2). To show that it also satisfies the condition Eq. (1.2) 
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we perform the fohowing lengthy calculations: Let Y^, G ^{Q) be arbitrary horizontal vector 
fields, i.e., Y^, Z^ G Vg for any q £ Q. We start from the following identity: 

djiY^,Z^) = Y^jiZ^)] - Z^^bCy^)] - 7{[Y'',Z^]). (4.4) 

The goal is to show that the right-hand side vanishes. Let us first evaluate the first two terms on 
the right-hand side of the above identity at an arbitrary point q £ Q: Let Zq := TqTTQ{Z^) G TqQ, 
then = hl^(Zg). Thus, using Lemma A. 2l we have 



f{<i) 

Hence, defining a function 72 : Q — ^ by 



liZ'^m = {\i\^omyfodW{q)M^{Zq 
= {^yfodW{q),Zq) 

^ dW{Z){q). 



lz{q) ■.= ^^dW{Z){q), 



we have '^{Z^) = 7^ o vr. Therefore, defining Yq := TqiiiY^^), i.e., Y^ = hl^(l^), 
Y%{Z'^)]{q)=Y\^zo^]{q) 

= (d(7zo7r)g,y,^) 



(7r*d7z)„n'' 

d^z{q),Tq^{Y^) 
{d^ziq),Yq) 

ybzm 

"1 



Y 



^ dW{Z) 



jY[Z[W]]-j^dfiY)dW{Z))iq). 



Hence we have 



Y^[^{Z^)] - Z^[^{Y^)] = y{Y[Z[W]] -Z[Y[W]])- j^{df{Y) dW{Z) - df{Z) dW{Y)) 



dW{[Y,Z]) - ^df AdW(Y,Z), (4.5) 



= ^^dWi[Y,Z])-j 

where we have omitted q and q for simplicity. 

Now, let us evaluate the last term on the right-hand side of Eq. (4.4): First we would like to 
decompose [y'^,Z'^]g into the horizontal and vertical part. Since both Y^ and Z^ are horizontal, 
we havtEl 

hori[Y\Z%) =hlf {[Y,Z]q), 

whereas the vertical part is 

ver([y^ ZX) = (M[Y\ zX)) (q) = - (s,(yi^, Z^)) ^ {q), 



^Recall that / : T*Q — >■ R is fiber-wise constant and thus, with a slight abuse of notation, we may write f{aq) = 
f(q) for any a, £ T^Q', therefore / may be seen as a function on Q as well. 
^See, e.g., Kobayashi and Nomizu [T71 Proposition 1.3 (3), p. 65]. 
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where we used the following relation between the connection A and its curvature B that hold for 
horizontal vector fields and Z^: 

= Y^iAiz'^Mq) - - A{[Y\ z'^m 

= -Ai[Y\z'^mq). 
As a result, we have the decomposition 

[Y\ z\ = hi^([y, Z],) - {b,{y^, Z,^))^ {q). 

Therefore, 

7([r^Zh])(<?) = {U^o^>,/fodWo^{q)M'^{[Y,Z],)) 

hlf oM/i/^ o dW o 7r{q),(B,{Y^, Z^)) ^ {q) 



^1// o dW{q), [y, Z]q) - (j(hlf oxl/i/^ o dW{q)) ,B,{Y^, Z, 



-^^dW{[Y,Z]){q) - (Johlf (dT^(g-)) ,S,(hl^(y5),hl^(Z,))> 



d7(y^z^) = —-^df^dW{Y,z) + -{dW)*E{Y,z) 



= ^^dW{[Y,Z]){q) - -j^^{dWrE{Y,Z){q), (4.6) 

where the second equality follows from Lemma |A.2| and the definition of the momentum map J; 
the fourth one follows from the linearity of hi and also of J in the fiber variables; the last one 
follows from the definition of H in Eq. (2.12): Since ttq o dW = idg and thus Tttq o TdW = idrpq, 
we have 

{dWrE{Y,Z){q) = E,^^^^{TdW{Y^),TdW{Zg)) 

= {Johl^{dWiq)) ,B,{hlfiY^,hlfiZ,))) . 
Substituting Eqs. (4.5) and ( |4.6[ ) into Eq. (4.4), we obtain 

j^df AdW(Y,Z) + j, 

= -j^{dfAdW-f (dWyE) (Y, Z) 

= -j^{dWr{dfAQ-fE) {Y,Z) 
= 0, 

where the third line follows sinc(|^(dT^)*/(g) = f[dW{q)) = f{q) and also that {dW)*@ = dW [see, 
e.g.,[Tl Proposition 3.2.11 on p. 179]; the last line follows from Eq. (3.10), which is assumed to be 
satisfied. □ 

5. Examples 

Example 5.1 (The vertical rolling disk; see, e.g., Bloch Consider the motion of the vertical 
rolling disk of radius R shown in Fig. [2j The configuration space is 

Q = SE{2) X §1 = {SO{2) K M^) x = {{ip,x,y,ij)}. 

Again recall that / : T'Q — > K may be seen as a function on Q as well. 
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Figure 2. Vertical rolling disk. 



Suppose that m is the mass of the disk, / is the moment of inertia of the disk about the axis 
perpendicular to the plane of the disk, and J is the moment of inertia about an axis in the plane 
of the disk (both axes passing through the disk's center). The Lagrangian L : TQ — )• M and the 
Hamiltonian H : T*Q are given by 

1 /.O .0\ It.'? 1 

2 



L = ^m[x^ + y2) + - J(^2 ^ 



and 



^=2^^ + T + T 



m J 
The velocity constraints are 

X = R cos ifijj, y = R sin if tjj, 

or in terms of constraint one-forms, 

u}^ = dx — R cos (f dil), uP' = dy — R sin 93 dil^. 

So the constraint distribution T) C TQ and the constrained momentum space M C T*Q are given 
by 

V = I ((^, X, y, tp) G TQ \ X = R cos -0, y = R sin V"! 

and 

mi? mR 



■M. = |(p<^,Px,Pt/,p^) e T*Q I p^. 



COS(/?P^, Py 



Let G 



and consider the action of G on Q defined by 

G xQ ^ Q; ((a, b), {cp, x, y, -0)) ^ {ip,x + a,y + b, ip). 



Then, the system is a Chaplygin system in the sense of Definition 2.1 The Lie algebra g is identified 
with M? in this case; let us use (^, ?/) as the coordinates for g. Then, we may write the connection 
^ : TQ g as 

d d 
A = (dx — R cos (p dtp) fS) 77- + (dy — R simp dtp) (8) t— , (5-1) 

at, or] 

and hence its curvature as 



sin ip dip A dip — — cos ipdip A dip ^ — 
ot, or] 



(5.2) 
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Furthermore, the momentum map J : T*Q ^ q* is given by 

J{Pq) =Px dC+Pydr]. 
The quotient space is Q := Q/G = {((/?, V')}- The reduced Hamiltonian H : T*Q — >■ M 
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IS 



H 



1 /I 



A simple calculation shows that the horizontal lift hi : T*Q — t- 7W is given by 



M . 



hl'^(p</p,Pv) = ( Pv^ 



mR 



mR 



(5.3) 



(5.4) 



(5.5) 



Then, we find from Eq. (2.12) along with Eqs. (5.1), (5.2), (5.3), and (5.5) that H = 0. Therefore, 
the sufficient condition, Eq. (3.10), for Chaplygin Hamiltonization reduces to d/A0 = 0, and hence 
we may choose / = 1. Thus, the Chaplygin Hamiltonian Hq : T*Q — t- M is identical to H (see 
Eq. (pf). 



To illustrate Theorem 4.1, we begin with the Chaplygin Hamilton-Jacobi equation (4.1): 

2' 



IfdWy I + mRWdW\ 



J \ dip 



+ 



P 



dip J 



E. 



(5.6) 



Now, we employ the conventional approach of separation of variables, i.e., assume that W : Q 
takes the following form: 



Then, Eq. (5.6) becomes 



IfdW^V I + mR'^ / dW, ^ ^ 



(- 

J \ dip 



+ 



P 



dip 



E. 



Since the first term on the left-hand side depends only on ip and the second only on ip, we obtain 
the solution 



dW^ 
dip 



7° 



dW, 



dtp 



1% 



(5.7) 



where 7^ and 7^ are the constants determined by the initial condition such that 



E. 



Then, Eq. (4.3) gives 



l{<f,x,y,tp) =-f^dip + 



mR 



cos 



ap^^^dx ^ — — sm 



'Pl4,dy + 7'^dip, 



(51 



which is the solution of the nonholonomic Hamilton-Jacobi equation ( 1.1 ) obtained in Ohsawa and 
Bloch [251 Example 4.1]: 

Example 5.2 (The knife edge; see, e.g., Bloch |3]). Consider a plane slanted at an angle a from 
the horizontal and let (x, y) represent the position of the point of contact of the knife edge with 
respect to a fixed Cartesian coordinate system on the plane (see Fig. [3]) and ip the angle of it as 
shown in Fig. |3j The configuration space is 

Q = SE{2) = 50(2) K = {i^,x,y)}. 
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Figure 3. Knife edge on inclined plane. 

Suppose that the mass of the knife edge is m, and the moment of inertia about the axis perpendicular 
to the inclined plane through its contact point is J. The Lagrangian L : TQ — )• M and the 
Hamiltonian H : T*Q are given by 

L = -m(i:^ + y^) + + mgxsina 

and 

„ i(pI + pI pI\ 

H = - \ \ — ;- — mgx sm a. 

2\ m J I 



The velocity constraint is 

and so the constraint one-form is 



sin ipx — cos (fy = 



u!^ = sin ipdx — cos 99 dy. 

The constraint distribution D C TQ and the constrained momentum space A4 C T*Q are given by 

V = {{if, x,y) G TQ I sin ipx — cos ipy = 0} 

and 

■M. = {{p^,Px,Py) G T*Q I sin(^p^ = coscppy} . 
Let G = M and consider the action of G on Q defined by 

GxQ^Q; {a,{(p,x,y)) ^ {ip,x,y + a). 



Then, the system is a Chaplygin system in the sense of Definition [2Tj The Lie algebra g is identified 
with M in this case; let us use 77 as the coordinate for g. Then, we may write the connection 
A: TQ ^ Q as 

d 

A = (dy — tanip dx) —, (5.9) 

or] 

and hence its curvature as 

1 d 

B= 7^dxAdip(^—, (5.10) 

cos^ if or] 

where we assume that ip stays in the range (0, 7r/2) or (7r/2,7r) to avoid singularities. Furthermore, 
the momentum map i :T*Q ^ q* is given by 

^{Pq)=Pydr]. (5.11) 
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The quotient space is Q := Q/G = {(99, x)}. The reduced Hamiltonian H : T*Q 



IS 



H 



1 / cos^ ip 



mgxsma. 



2 \ m 

A simple calculation shows that the horizontal lift hl''^ : T*Q ^ Ai is given by 

hl^ {p^,Px) ={p^,cos'^ ippx, simp cos ip px) ■ 



Then, we find from Eq. ( |2.12[ ) along with Eqs. (jS^, ( |5.10[ ), ( |5.1l[ ), and ( |5.12[ ) that 

S = tan ipdx A dip. 



Therefore, the sufficient condition, Eq. (3.10), for Chaplygin Hamiltonization gives 

dl 

' dx 



It is easy to find the solution 

/ = cos ip. 

Note that / is nowhere-vanishing if ip is assumed to be in the range (0,7r/2) or {tt/2,tt). 
Then, Eq. (3.8) gives the following Chaplygin Hamiltonian: 

Pv> Px 



Hc{ip,x,p^,Px) = H[ip,x 



cos if COS if 



1 / 1 



pi + 



1 2 

J COS^ ip ^ 



mgx sma. 



The Chaplygin Hamilton-Jacobi equation (4.1) then becomes 

2 



1 fdWV 
dx J 



+ 



1 



fdW 

\dip 



Assume that W : Q 



J cos^ ip 

takes the following form: 

Wiip,x) = W^iip) + Wxix). 



mgxsma = E. 



Then, Eq. (5.14) becomes 



1 fdW, 



m\ dx 



{2mgsm.a) x + 



1 



dW^ 



J cos^ ip\ dip 



E. 



(5.12) 



(5.13) 



(5.14) 



The first two terms in the brackets depend only on x, whereas the third only on 99, and thus 



1 fdW,. 



m\dx 



{2mg sin a)x = 2E 



1 



J ' 



dW^ 



cos^ if \ dip 



0^2 
ip) 1 



with some positive constant 7^. Hence, assuming dWx/dx > 0, we have 



dWx 
dx 



L(2^-<f! 



+ [2m'^ g sm a) X , — 

dip 



1^ cos ip. 



Then, Eq. (4.3) gives 



7((/.,x,y)=7°d(/.+ ^/m(2i?-^ 



+ {2m'^g sin a) x (cos ipdx + sin ipdy), 



which is the solution of the nonholonomic Hamilton-Jacobi equation ( 1.1 ) obtained in Ohsawa and 
Bloch [251 Example 4.2]. 
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6. Further Reduction and Hamiltonization 
It often happens that there does not exist an / that satisfies the necessary and sufficient condition 



in Proposition 3.9 or the sufficient condition, Eq. (3.10), and hence we cannot Hamiltonize the 
system based on the above theory. However, we may reduce such systems further and then attempt 
to Hamiltonize the further-reduced system. 

6.1. Further Reduction of Chaplygin Systems. We consider the fohowing special case of the 
"truncation" of Hochgerner and Garcia-Naranjo [15] Section 3.B]. Recall the reduced Chaplygin 



system, Eq. (2.10), on T*Q, i.e., 

ixO.''^ = dH (6.1) 

with the almost symplectic form 

f^nh .= f7 _ (6.2) 

and consider a free and proper Lie group action K x Q ^ Q, oi : Q ^ Q with any k ^ K, that 
satisfies the following conditions: 

I. The Hamiltonian H is i^-invariant, i.e., H o T*^^j^ = H for any k ^ K, where T*^^ is the 
cotangent lift of 

II. For any element r] in the Lie algebra t of K, the infinitesimal generator i]j<,q satisfies 

VqS = 0- (6-3) 

Now, let Jk : T*Q ^ i* he the equivariant momentum map for the cotangent lift of the i^-action 
i.e., for any aq ^ T*Q and rj & t, 

(Jx(ag),r?) = (ag,r/Q) • (6.4) 
Also define : T*Q — )• M by J^{aq) := {J Kictq),v) each tj G t. Then, we have 



Notice that Condition [Tl] implies 
and thus 



„• o- dp 



In other words, Jk is a momentum map with respect to both the standard symplectic form Q and 
the almost symplectic form Cl""^. We also have the following: 



Proposition 6.1. Under Conditions'^ and\l^ stated above, the momentum map Jk '■ T*Q — )• t 



IS 



conserved along the flow of the vector field X of the reduced Chaplygin system, Eq. (6.1). 
Proof. Follows easily from the following calculation: 

X[Jl] = ^^dJl 

= -ir,^,^ix^'^^ 

= -hT*QdH 

= -^T'q[H] 
= 0, 

where we used Eq. (6.5) in the second line, and Condition |l] in the last line. □ 



Also, let Kfj^ be the coadjoint isotropy group of i.e., := {k G K \ Ad^ ^ = /i}, and assume 
III. /Lt G F is a regular value of Jk, and acts freely and properly on J^{^). 
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Since 3k is a momentum map with respect to the ahnost symplectic form Q"^, the two-form 
0°*^ itself works as a "truncated form" (see Hochgerner and Garcia-Naranjo [15', Section 3.B and 
Theorem 3.3]) in this special case: Performing the almost symplectic reduction of Planas-Bielsa 
[27], we may drop the dynamics to 3j}{fi)/Kf^ as follows: 

Proposition 6.2 (Further Reduction of Chaplygin Systems). Under Conditions^ Ilh we have the 
following: 

(i) There exists an almost symplectic form Vl"^ on 3j}{fi)/K^ uniquely characterized by 

7rX' = i;^^°', (6-6) 
where ifj, : J^^(^) ^T*Q and vr^ : J^^(/i) 3j}{fj,)/ Kfj_. 
(a) The reduced Chaplygin system, Eq. (|6.1|), is further reduced to the following system: 



i^Qf = dH,, (6.7) 
where X and are n ^-related, i.e., 

Ttt^ o X = X^o TTf,, (6.8) 
and Hfj_ : 3^{ii)/K^ — )• M is defined by 

H^o7r^ = H oi^. (6.9) 



(Hi) The almost symplectic form is written as 
where H„ is uniquely characterized by 



nnh _ Q _ - 



7r;E, = i;E. (6.10) 



Proof. ^ and ^ follow directly from Planas-Bielsa [27^ Theorem 2.1]. (|in]) Since 3k is an equi- 
variant momentum map with respect to the canonical symplectic form Cl, the symplectic reduction 
of Marsden and Weinstein [23] applies here as well (not to the reduction of the dynamics but to 
the reduction of the symplectic structure). Hence there exists a unique (strictly) symplectic form 
Clfj^ on 3^{n)/Kfj, such that 

vr*^^ = il^. (6.11) 

Combining this with Eq. (6.6), we have 



Since vr^ is a surjective submersion, the pull-back vr* is injective, and thus the uniqueness follows. □ 

Furthermore, under certain assumptions, we may employ a result from the theory of cotangent 
bundle reduction (see, e.g., Marsden et al. Section 2.2]) to make our result more explicit. To 
that end, we first define a mechanical connection on the principal bundle 

7f : Q ^ Q/K =: Q 

as follows: For each q £ Q, let I{q) : € — t- £* be the locked inertia tensor defined by 

{l{q)v,0 = 9q{'nQ{q),QQ{q)) , 

where g is the kinetic energy metric defined in Eq. ( |2.7[ ), and r/ and C, are arbitrary elements in t. 
Then, the mechanical connection Ak ■ TQ — t- t is defined by 

AKivg) ■■= liqy^ ° 3K{¥L{vii)) . (6.12) 
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We will also need the "/x-component" of Ak, i-e., the one- form on Q defined by a^(g) := 
^if(g)*/U, or equivalently, 



Let us introduce the two-form (3^ on Q defined by 



and also the two- form on T*Q defined by 



B 



K ._ 



(6.13) 
(6.14) 

(6.15) 



where tTq : T*Q ^ Q is the cotangent bundle projection. 
Now, we assume the following: 

IV. Kfj^ = K, which is always the case when K is Abelian; 
V. is X-invariant and takes values in J^^(^). 



With these additional assumptions, we have the following important special case of Proposition 6.2 



Proposition 6.3 (Further Reduction of Chaplygin Systems — Special Case). If, in addition, Con- 
ditions IV_ and\^hold, then we may extend the results of Propositioned^ so that the dynamics after 
the second reduction is described on the cotangent bundle T*{Q/K) = T*Q as follows: 

(i) The reduced space 3j}{fi)/K is symplectically diffeomorphic to T*Q with the symplectic struc- 
ture O — , where O is the standard symplectic form on T*Q. 

(a) Let (ffj_ : 3^{pL)/K T*Q be the symplectomorphism that gives the correspondence in 
Then, the dynamics on 3~f^{^)/K defined by Eq. (6.7) is equivalent to the one defined by 



nh 



iv ^ 



dH., 



where 



with X 



n-^:=n-Bj^ 



(6.16) 
(6.17) 



M — {^f.^)*^t^> % iv^^y^n., and H, 



Proof. ^ See Marsden et al. [24t Theorem 2.2.3 on p. 64]. The construction of the map 99^ is 
summarized in Appendix |B| The diagrams below summarize the spaces and almost symplectic 
forms involved in the procedure of the reduction. 



T*Q 

A 



{l^)/K 



T*Q 



Of h 



^ Apply (v?^^)* to both sides of Eq. ([6^ and use the fact from ^ that {^p-^)*nf, = O - Bjf . □ 
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6.2. Hamiltonization after Second Reduction. Now, we follow a similar argument as in Sec- 
tionjsjto discuss the Hamiltonizability of the system defined by Eq. (6.16): Let : T*Q — )• M be a 
smooth nowhere- vanishing function that is constant on each fiber, and define the map : T*Q ^ 
T*Q by 

% : « ^ //^ «• 



Define the vector field analogously to Eq. (3.2) so that 



and hence X^//^ and Xj^ are ^'/^-related: 



(6.18) 



(6.19) 



Following the same arguments as in the proofs of Proposition 3.9 and Theorem 3.13 we obtain 
similar results for the further-reduced system, Eq. (6.16). 

Proposition 6.4 (Necessary and Sufficient Condition for Hamiltonization after Second Reduction). 

The vector field X^ G X(T*Q) satisfies Hamilton's equations 



1 -CrU,^ 



with the Chaplygin Hamiltonian 
if and only if the one-form i 



df,Ae-f^jBff + ^ E, 



(6.20) 

vanishes, where Q is the symplectic 



one-form on T*Q. 



Proof. The result follows from essentially the same calculations as in the proof of Lemma 3.8 The 
only difference is the treatment of the curvature term Bjf , which is not present in Lemma 
Specifically, we need to calculate Bjf: From the definition of Bjf , Eq. (6.15), we have 



3.8 



IT/ 



Q ^ 



where we used the fact that ^i//^ is fiber-preserving, i.e., tTq o '^i/f^ 

1 r 



TT/s- Therefore, we obtain 



n + 



df,Ae-f'^(Bj!+^l^fE, 



and thus the claim follows. 



(6.21) 



□ 



Remark 6.5. Since Q — Bjf is also a (non-standard) symplectic form as well, we may discuss 
Hamiltonization with respect to this symplectic form. However, we prefer to work with the standard 
symplectic form since the standard Hamilton-Jacobi theory directly applies to Hamiltonian 
systems defined with the standard symplectic form 0,. 

Theorem 6.6 (A Sufficient Condition for Hamiltonization after Second Reduction). Suppose there 
exists a nowhere-vanishing fiber-wise constant function f^ : T*Q — t- M that satisfies the equation 



df,Ae = f'^(Bj^+^l/^E, 



(6.22) 
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Then, the vector field € X(T*Q) (see Eq. (6.18 )j satisfies the following Hamilton's equations: 



(6.23) 



and, as a result, the further-reduced nonholonomic dynamics, Eq. (6.16), has the invariant measure 
/"~^A, where n := dimQ and 



A := ^ J7 A • • • A O . 



n! 



Proof. Follows immediately from Eq. (6.21) and Corollary 3.4 



□ 



7. NONHOLONOMIC H-J THEORY VIA HAMILTONIZATION AFTER SECOND REDUCTION 

7.1. Relationship between the Chaplygin H J and Nonholonomic H J Equations after 



Second Reduction. If the system is Hamiltonized in the sense of Theorem 6.6 then we have the 
Chaplygin Hamilton-Jacobi equation 



Hi^ o dWi" 



E 



(7.1) 



corresponding to Hamilton's equation (6.23). One then wonders if there is any relationship between 



4.1 



W^^ and 7 that is similar to the one obtained in Theorem 

A natural starting point towards the answer to this question is, again, to look into the relationship 
between the Chaplygin Hamiltonian and the original Hamiltonian H; then we obtain the 
relationship between the two solutions and 7 by exploiting the geometry involved in the process 
of reduction and Hamiltonization. The diagram below combines the following things together: the 
first and second reductions of Chaplygin systems; the relationship between the two Hamiltonians 



J^^(^) — > J^^(O) (see Appendix B|; also the horizontal lift 



Hq and H] the shift map shifty 

hl-^ ■.T*Q^M := J ^^(O), which is defined in a similar way as hl^ (see Eq. ^ ) using the 
connection Ak (see Eq. (6.12)) as follows: Let us define the horizontal space 



V := ker^ii- C TQ. 



Then, the connection Ak induces the horizontal lift hl^ : TQ 
Let us also define 



V defined by hi 



V 



(Tvfl 



V) 



M:=J^i(0) 



Then, it is straightforward to see that M. = FL(D). Now, we define the horizontal lift hl-^ : T*Q 
Ai as follows: 



hlf :=FZ,-ohf o(FZ)ri, 



(7.2) 
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where L:TQ is defined by L := L o hP. 




(7.3) 



T*Q 



That the map hl^^ fits into the diagram is shown in Appendix|c| (see also Appendix b|. The diagram 
also shows the map dW^ : Q ^ T*Q with being a solution of the Chaplygin Hamilton-Jacobi 



equation (7.1); this leads us to the following result that is similar to Theorem 4.1 



Theorem 7.1. Suppose that there exists a nowhere-vanishing fiber-wise constant function : 
T*Q —7- M that satisfies Eq. (6.22), and hence by Theorem 6.6, we have Hamilton's equations 



( |633l ) for the vector field X^. Let W : Q 
equation (7.1), and define 7 : Q — t- 7W by 



R be a solution of the Chaplygin Hamilton-Jacobi 
7(g):=hl^o7^ovr(g) (7.4) 



with j/j, : Q ^ T*Q defined b 



shifty ^ o hl^ o'^i/f^ o dW o Tt{q) 



io ° hi; 



M 



ffj, 



dW''{q)]+a^{q), 



(7.5) 



where q := 7r(g) and q := vr(g). Then 7 satisfies the nonholonomic Hamilton-Jacobi equation (1.1) 



as well as the condition Eq. (1.2). 



Proof. That the one- form 7 defined by Eqs. (7.4) and (7.5) satisfies the nonholonomic Hamilton 



Jacobi equation (1.1) follows from the diagram (7.3). Showing that it also satisfies the condition 
Eq. (1.2) requires tedious calculations: Following a similar calculation to that of 1^7 (1"^, Z^) in the 



proof of Theorem 4.1, Eq. (7.4) gives 

dj{Y'^,Z'^) = d%{Y,Z) + %^{Y,Z) 
for arbitrary horizontal vector fields Y^, G ^{Q), where Y := T'ir{Y^) and similarly for Z. 



(7.6) 



Let us calculate the first term in Eq. (|7.6|): Writing 

1 



70 := io o hi 



M 



dW^ 



^Qee Appendix [b| for the relationship between i,,, io, and shifty: We have o shift,j^(pq) — io{Pq) + citiil) for any 

p.-e (0). 
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we have 7^ = 70 + and thus 

d%{Y, Z) = d%{Y, Z) + da^{Y, Z). 

Calculation of d^y^iY^Z) is somewhat similar to that of d'^iY^ , Z^) in the proof of Theorem 
but there is one difference: Y and Z are not horizontal here. Specifically, we have 

d^o{Y, Z) = y[7o(^)] - Z[^^{Y)] - 7o([y, Z]) 



4.1 



1 

fu 



1 



dW^{[Y, Z])--^df^A dW^iY, Z) - 7o([y, Z]), 
J I 



where we defined Y := T7r(Y) and similarly for Z. To calculate 7o([y, .Z']), we decompose [Y, Z] 
into the horizontal and vertical parts: 

[Y,Z] = hf{[Y,Z])+iAKi[Y,Z]))Q, 

where we note that TTt{\Y, Z]) = [Y, Z], since Y and Z are vf-related to Y and Z, respectively. As 
a result, we have 



7o([y,z])(9) 



1 



dWMy,Z]){q) + 
dW^{[Y,Z]){q), 



1 



Jk o hlf ( dW^ {q)),AK{ [Y, Z] ) (q) 



where the second term vanishes because hi takes values in := J^, (0). Next let us calculate 
da^{Y,Z): Using Eq. (6.14), the relation tTq o dW^ = idg, and Eq. (6.15), we obtain 

da^{Y,Z) = Tt*/3^{Y,Z) 

= r o{7TQodWn*(3^.{Y,Z). 

= r o(dTy^)*ovrt/3,.(y,z). 

= Tt* o{dW^'yBff{Y,Z). 
= {dWf')*Bl^{Y,Z). 



Therefore, the first term on the right-hand side of Eq. (7.6) becomes 

d%{Y, Z) = -^{dWn* (df, A e - /2 Bj^) {Y, Z), 



since (dW^)*/^ = /m and also {dW^,)*e = dWf,. 

Now, let us evaluate the second term on the right-hand side of Eq. ( 7.6 ) : Substitution of Eq. ( 7.5 ) 
gives 

-1 „ u^M 



7u" 



shifty' o hlf" ol-i/^^ o dW o 7f 
shifty 1 o hlf o§i/ o dW^ o Tf) ' o tt;-, 
ir^ o shift^^ o hlf o^yj^ o dW^" ottYe 



1//m 



dW^ 



o vr 



vr,, o shift,, ^ o hl^ 



^i/f^ o dW^^ o vr 



O dWf" o 7f 



7:*o{dWn*o^l/^E, 
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where we used Eq. (6.10), the relation vr^ o shifty ^ o hl^ 
definition of H„ from Proposition 6.3 This imphes 



from the diagram (7.3), and the 



%E{Y, Z) = {dw^^r o §1/. Z). 



As a result, Eq. (7.6) becomes 



1 

72' 
J ^ 



which vanishes because the sufficient condition, Eq. (6.22), is assumed to be satisfied. 



□ 



7.2. Example of Further Reduction, Hamiltonization, and Chaplygin H J Equation. 

Example 7.2 (The Snakeboard; see, e.g., Ostrowski et al. [26j, Bloch et al. [1] and Koon and 
Marsden [20J ) . Consider the motion of the snakeboard shown in Fig. |4j Let m be the total mass 




Figure 4. The Snakeboard. 

of the board, J the inertia of the board, Jq the inertia of the rotor, Ji the inertia of each of the 
wheels, and assume the relation J + Jq + 2 Ji = mr^. The configuration space is 

Q = SE{2) X §1 X §1 = {SO{2) X M^) x x = {{e,x,y,(f),ij)}. 

The Lagrangian L : TQ — t- M and the Hamiltonian H : T*Q are given by 

1 



L 



2 L 



m[±^ + y' + r'e'') + 2Jo V + 2 Ji + Jo V'' 



„2fl2 



and 



^ 2m^^^'^^^^ ^ 2(mr2- Jo) 



The velocity constraints are 

X + r cot (j) cos 9 9 = 0, y + r cot cj) sin 99 = 0, 

or in terms of constraint one-forms, 

u!^ = dx + r cot <p cos dO, uP' = dy + r cot (p sin 9 d9. 

So the constraint distribution T) C TQ and the constrained momentum space Ai <ZT*Q are given 
by 

P = G TQ I u:\9,±,y,^,i;) = 0,s = l,2Y 

and 

M = {{pe,Px,Py,P4>,P'^) G T*Q I = -Kcotc/) cos9{p0 -p^), py = -Kcot(j) sm9{p0 - p^p)} , 

where k := mr/{mr^ — Jq). 

Let G = 'M? and consider the action of G on Q defined by 

GxQ^Q; ((a, b), {9, x, y, (p, ^)) ^ {9,x + a,y + b, (p, ip). 
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Then, the system is a Chaplygin system in the sense of Definition 2.1 The Lie algebra g is 
identified with in this case. Let us use again t]) as the coordinates for g. Then, we may write 
the connection A : TQ — )• g as 



d d 
A = (dx + r cot 6 cos 6 dO) 1X177- + idy + r cot 6 sin dO) ® — , 

ot, or] 



and hence its curvature as 
B 



2 d o d 

r cos 6 CSC (p dO A d<j) (S> + r sin 6 esc (pd9 Ad(p <^ — 

at, or] ^ 

Furthermore, the momentum map J : T*Q — )• g* is given by 

J(Pg) =Px dC+Pydr]. 

The quotient space is Q := Q/G = {{9, (j), ip)}, and the reduced Hamiltonian H : T*Q 



(7.7) 

(7.8) 

(7.9) 

: is 



H 



sm 



Jn sin^ 1 



{Pe - Pv) + 



2 , P<l> , 



2Ji Jo 



and the horizontal lift hl^ : T*Q ^ A4 is given by 



(mr^ — Jo) sin^ < 
mr'^ — Jo sin^ d 



-ipe - Pi,) 



mr cos (p sm (p cos ( 
mr^ — Jo sin^ (p 
mr cos d> sin 6 sin 1" 



(Pe-P^), 



mr^ 



Jq sin 



Then, we find from Eq. ( [2l2| ) along with Eqs. ((T^T]), ([T^, ([T^, and ((Tlo]) that 

„ mr^ (pe — p^ ) cot 



mr 



2 - Jo sin^ I 



^6* A del). 



However, there exists no function / that satisfies the sufficient condition, Eq. (3.10), for Chaplygin 
Hamiltonization. In fact, one can show (see p^) that there does not exist an / which satisfies 
the necessary and sufficient condition for Hamiltonization from Proposition 3.9 Hence the system 
is not Hamiltonizable at this level of reduction. Therefore, we would like to further reduce the 
system: Let K = 'B^ and consider the action of -ftT on Q defined by 

KxQ^Q; {c,{9,^,tP))^{0,ct>,i; + cy, 

and so ^^{9, 4>, ip) = (9, (j),ip + c) for any c ^ K. This gives rise to the cotangent lift 

K X T*Q ^T*Q; {c, {9, (I),il),p0,p^,p^)) ^ {9, 4>,ip + c,pe,p^,p^), 

that is, 

T*^^^{9,cl),'tp,pe,p^,p^) = {9,(/),'ip + c,pe,p^,p^). 
It is easy to see that the Hamiltonian H is ii'-invariant. Also, for any G t,we have the infinitesimal 



generator Ct*q — Cd/dip and so easily see that i^^ 



0. Hence Conditions 



Therefore, by Proposition |6.1[ the corresponding momentum map 

Ji^(Pg) =Pi^dC 



|T|and|H| 



are satisfied. 



is conserved. It is straightforward to check that Condition |III| is satisfied for any n = dC, G V . 
Then, Proposition 6.2 gives the reduced dynamics on 3^{pL) / K^, and Eqs. (6.9) and (6.10) give 



mr^ 



. ■ 2 +VT^T 

— Jq sm (p IJ\ Jo 
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and 



„ _ mr'^ipg - fi^) cot(j) 
^ mr^ — Jq sin^ (p 

Furthermore, Eq. (6.12) gives the mechanical connection 

Ak = [dO + dilj) (g) 



dO A dcj). 



d_ 



(7.11) 



and hence Eq. (6.13) gives 

a^ = fi{dd + dij), (7.12) 
and so = and Bjf = 0. It is also straightforward to check that Conditions IV and|v]are satisfied. 
Therefore, we may apply Proposition 6.3 to this case. Specifically, we have Q := Q/K = {{9,(j))}, 
and Eq. ( B.5[ ) (from Example B.l in Appendix [B]) gives 

ipj,^ : T*Q J]^^{n)/K; {e,(j),pe,p<i,) H> {e,(l),pe + 

and hence we have 



-,-1/ 



and 



H^{e,4>,pe,p^) := Hf, o (f^-[(),(p,pg, 



sm 



1 



2 \ mr"^ — Jo sin 



2 , , 



2Ji ^ J( 



mr pg cot (j) 
mr'^ — Jq sin^ 



d9 A #. 



Therefore, the sufficient condition, Eq. (6.22), for Chaplygin Hamiltonization becomes 

mr^ cot (j) 



-P9 



which gives 



de 



dfa 



A straightforward integration yield^ 



mr'^ — Jq sin^ ( 
mr^ cot (j) 

~ 9 T ■ 2 1 Jf^' 

mr'^ — Jq sm (p 



\l mr"^ — Jq siv? (f) 



where we assume that | sin(f>\ < ^mj Jqt. Then, Eq. (6.20) gives the following Chaplygin Hamil- 
tonian: 



H^{e,^,pe, 



\/ mr"^ — Jq siT? (j) a/ mr^ — Jq sir? < 



2 7 c;v.2 



1 / 2 ~ Jo sin 



2 Ji sin^ 



■Pe, 



sm ( 



Hence the Chaplygin Hamilton-Jacobi equation (7.1) becomes 



de 



+ 



mr — Jq sin 
2 Ji sin^ 4) 



^W^' 



+ 



/4 



E. 



(7.13) 



Assume that : Q — s- M takes the following form: 



^For mr^ = Jo = 1, this verifies the result of [HI Section 4.4]. 
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Then, Eq. (7.13) becomes 



+ 



mr 



2- Jnsin^c/, fdWjf' 



2 Ji sin^ 



+ 



f4 

Jo 



E. 



The first term in the brackets depends only on 9 whereas the second only on cp, and the third one 
is constant. Thus we have 



dW!f 



sm ( 



dcj) ^ jnr"^ - Jq sin^ i 
with some set of constants and 7^ that satisfy 

which is solved for 7^ (assumed to be positive) to give 



7^, 



E, 



\ 



2 E 



ill? ^\ ' 
4J1 2Jo 



Therefore, Eq. (7.4) with Eq. (7.5) gives 

(mr^ — Jq) C sin ( 



7(6*, X, y, ^, V') = /x^ + 



mrC cot (j) sin (p 



de 



(cos 6 dx + sin 6 dy) + 7^ dc/i + fi^ dip, 



where we defined 



C :-- 



ill? ^% 



'E 



4Ji 2 Jo ' 



9K 



:= \J {mr'^ — Jq sin^ (/)) /2. 



This is the solution of the nonholonomic Hamilton- Jacobi equation (1.1) obtained in Ohsawa and 
Bloch [251 Example 4.3]. 

8. Conclusion and Future Work 

We established a link between two different approaches towards nonholonomic Hamilton-Jacobi 
theory; the direct one in |16| |25] and the indirect one via Hamiltonization. We formulated the 
procedure of Hamiltonization in an intrinsic manner; this helped us understand the relationship 
between the two approaches and also lead us to the formulas relating the solutions of the two 
different types of Hamilton-Jacobi equations resulting from the direct and indirect approaches. 
The formulas provide us with the following new method to exactly integrate equations of motion 
of nonholonomic systems: 

1. Reduce and Hamiltonize the nonholonomic system. 

2. Solve the Chaplygin Hamilton-Jacobi equation for the Hamiltonized reduced system. 

3. Use the formula in Theorem 14.11 or 17.11 to obtain the solution of the nonholonomic Hamilton- 
Jacobi equation for the full dynamics. 

4. Integrate the full dynamics using the solution as shown in Ohsawa and Bloch [25j . 

A notable feature of this method is that it links the solution of the Hamilton-Jacobi equation for 
the reduced system with integration of the full dynamics. We illustrated this method with a few 
examples and obtained the solutions identical to those in Ohsawa and Bloch j25j . 
The following questions are interesting to consider for future work: 
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• Hamiltonization and Hamilton- J acohi theory for a more general class of nonholonomic sys- 
tems with symmetries: This paper only dealt with Chaplygin systems, a special case of the 
more general class of nonholonomic systems with symmetries treated in Bloch et al. . We 
are interested in extending our results to the general case, possibly relating them to the 
results on existence of an invariant measure in [30] . 

• Application to nonholonomic systems on Lie groups: Nonholonomic systems on Lie groups, 
such as the Suslov problem (see, e.g., Kozlov [21j, Zenkov and Bloch [29J), often involve 
an interesting question on integrability: Whether or not the full dynamics is integrable 
when the reduced dynamics is (see Fedorov et al. |13j). Relating this question with the 
Hamilton-Jacobi equations for the full and reduced dynamics is an interesting question to 
consider. 
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Lemma A.l. The horizontal lift hi is invariant under the action of the cotangent lift of 
Specifically, for any h £ G, we have 



Now, using the G-invariance of the Lagrangian L, we have, for any Vq G TqQ and w^q G T^qQ, 



Acknowledgments 



Appendix A. Some Lemmas on the Horizontal Lift hi- 




(A.l) 




{¥Lhq O Tq^h{Vq),Whq) 



d_ 
de 



L{Tq^h{Vq) +eWhq) 



£=0 



d_ 
de 



LoTq^h{Vq+eThq^h-l 



{Whq)) 



£=0 



= -^L{Vq + eThq^y,-i{Whq)) 

= {¥Lq{Vq),Thq^y,-l{Whq)) 
= {T;^h^.{¥Lq{Vq)),Whq), 



e=0 



and thus FL/^j o T^h = T*^fi-i o ¥Lq. Hence we obtain 



hl^=r;cI>^-ioFL,ohl^o(FZ)^r 



= r*ci>^-iohi^ 



□ 
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Lemma A. 2. Let q be an arbitrary point in Q and q = ir{q) E Q. For any aq E T^Q and Vq E TqQ, 
the following identity holds: 

{h\^{aq)M^{vq)) = {aq,Vq) . 



Proof. Follows from the definitions of g and hl^ (see Eqs. (2.7) and (2.8), respectively): 



{hl^ {aq)M^{vq)) = o hl^ o{^)--\aq)M^, {vq) 
= g,{yPqo{g%\aq)M'',{v,)) 



9q{i9^)q^{aq),Vq 



9\°{9^)q^{aq),v^ 
{aq, Vq) . 



□ 



Appendix B. Construction of : 3j^{^)/K T*Q 
We briefly summarize the construction of the map 99^ : J^^[fi)/K — t- T*Q that appears in 



Proposition 6.3 following Marsden et al. |24^ Section 2.2]. First define (fo : J^^(O) — )• T*Q by 

{(pQ{pq),TqTt{vq)) = {pq,Vq) (B.l) 

for any pq E TqQ and Vq E TqQ. Let ttq : J^^(O) — )• J^^(0)/i^ be the projection to the quotient. 
Then, ipQ : J^^(0)/i^ — T*Q is uniquely characterized by the relation 

(B.2) 

Proof of Theorem 2.2.2 



yjQ O VTo = LpQ. 



It can be shown that is in fact a diffeomorphism (see Marsden et al. 
on pp. 62-63]). We also introduce the shift map 



defined by 



Shifty ■.T*Q^T*Q 
Shifty (pg) := Pq - af,{q), 



(B.3) 



where is the one- form on Q defined in Eq. (6.13). This gives rise to the iC-equivariant diffeo- 
morphism 

shifty :J^i(/i) ^J^i(O), 
and the commutative diagram below, where in and io are both inclusions. 



T*Q 



Ja-(/") 



Shifty 



T*Q 



shift,, 



J7/(o) 



Since the map shifty is -fC-equivariant, it induces the diffeomorphism 

shifty : 3],Hl2)/K ^ 3kH0)/K. 
The map 93^ : J~^^{fi)/K — )• T*Q is then defined by 



fti-=fo° shifty 



(B.4) 
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The diagram below summarizes the construction of (p^. 

shift n 



■Jx(o) 



^0 



shift u 



T*Q 



Example B.l (The Snakeboard; see Example 7.2). Let us first determine ipo and (po. Note that 
we may parametrize 3j}{0) as follows: 



and also that 7r(^,0, ^) = {9,(t>) and hence TT:{ve,v^,v^) = {ve,v^). Therefore, Eq. (B.l) gives 

(fo{9,(l),ip,pe,p^) = {9,(l),pg,p^). 



Since 7ro{6,(t),ip,pg,p^) = {e,(f),pg,p^), Eq. ( |B.2[ ) gives 

ipQ{e,cl),pe,p^) = {e,(p,p0,p^). 



Now, let us determine the map shift^: Using the in Eq. (7.12), we find, from Eq. (B.3), 

Shih^{e,(j),ilj,p0,p^,p^) = {e,4>,ip,pe - fi^,p<i„p^ - IJ,^). 
Parameterizing J^^(/x) as 

(/^) = {{^,4>,i^,Pe,P4>,Pi,) G T*Q I p^ = 11^] = {{e,(f),i;,pe,p^)}, 

we obtain 

shiit^{e,(l),ip,pe,p^) = (6',(/),^/^,pe - fi^,p^), 

and hence 

shift f,{e,4>,p0,p^) = {9,(l),pe - fi^ 



As a result, we obtain, from Eq. (B.4), 



^fMiO,4>,Pe,P<t>) = (0,(1), pe - l^4„P<j,)- 

Appendix C. On the Horizontal Lift hl-'^ 

Lemma C.l. Let q be a point in Q and q = Tt{q). Then, we have hl^o(^o(Pg) 
Pq G 3^^{0) and also (pQ o hl^^ = idj,,Q, and the diagram 

hi-^ 



(B.5) 



Pq for any 



commutes with an appropriate choice of the base point q of the image of hl^ . 
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Proof. Let Pq ^ TiQ and Vq S TqQ be arbitrary. Then, Eq. ( |B.l ) imphes that 

(^Oohlf (p,-),r57f(hlf K-))) = {hlf{Pq),h\f{Vq 



'1' 

{Pq, Vq) , 



where we used an identit y on 
the same way as Lemma A. 2 



and thus the above equation 



pairings between the images of hl"^ and hl'^ , which can be shown in 



However, the definition hi := (T7r|-p) implies Ttt o hi = id^^g, 
reduces to 



(^ifO O hlf{pq),Vg^ = {pg, Vq) . 

Therefore, we have (po o hl"^ = id^.g. So Eq. ( B.2[ ) gives 

ipQO-KQO hl^ = id^.Q, 

which also implies 

^oohl^o(^o = idj-i(o)/i^, 

since (pQ is a diffeomorphism. 

To show hl^ °'fo{Pq) = Pq, take an arbitrary Vq G TqQ. Then, we may decompose Vq into the 
horizontal and vertical parts, i.e., 

Vq = hl?{Vq) + {AK{Vq))Q{q) 

where Vq = TqTt{vq). Therefore, we obtain 

|hlf 0^^{pq),Vq) = (hlf 0^o(P5), hlf (^,-)) + (hlf 0(^o(P5), {AK{Vq))Q{q)) 
= {<fo{Pq),Vq) + (Jk O hlf* O'P0{Pq),AK{Vq) 



= {<P0{Pq),TqTt{Vq)) 
= {Pq, Vq) , 

where we used the fact that hl"^ takes values in M ■= J]^^(0), and also Eq. (B.l). Hence 
hl^ o(fo{pq) = Pq and also, by Eq. ( |B.2[ ), hl^ oipo o -Ko{pq) = pg. □ 
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